Excitonic gap, phase transition, and quantum Hall effect in graphene 
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We suggest that physics underlying the recently observed removal of sublattice and spin degenera- 
cies in graphene in a strong magnetic field describes a phase transition connected with the generation 
of an excitonic gap. The experimental form of the Hall conductivity is reproduced and the main 
characteristics of the dynamics are described. Predictions of the behavior of the gap as a function 
of temperature and a gate voltage are made. 
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I. INTRODUCTION 

The properties of graphene, a single atomic layer of 
graphite have recently attracted a lot of attention, 
especially after the experimental discovery @,J1] and (in- 
dependently of that) theoretical prediction 0, Q of an 
anomalous quantization in the quantum Hall (QH) effect 
(for earlier considerations of the QH effect in graphene, 
see Ref. @). The graphene material is unique because 
of its band structure with two inequivalent Dirac points 
at the corners of the Brillouin zone. As a result, its low- 
energy excitations are described by effective "relativistic- 
like" Dirac equation where the speed of light is replaced 
by the Fermi velocity Vp Q. 

These relativistic-like features of graphene are at the 
heart of the anomalous integer QH effect. In this case, 
the filling factors are v — ±4(|n| + 1/2), where n is the 
Landau level index. For each QH state, a four-fold de- 
generacy takes place: it is the sublattice and spin degen- 
eracy for the lowest Landau level (LLL) with n = and 
the valley and spin degeneracy for higher Landau levels 
(LLs) with \n\ > 0. In the very recent experiments Q, it 
has been observed that in a strong enough magnetic field 
the new QH plateaus, v = 0, ±1 and ±4, occur, that was 
attributed to the magnetic field induced splitting of the 
LLL and the n = ±1 LLs. It is noticeable that while the 
degeneracy of the lowest LLL is completely lifted, only 
the spin degeneracy of the n = ±1 LLs is removed. 

In this paper, we suggest that the origin of the plateaus 
v = and, especially, v = ±1 is deeply connected with 
a phase transition with respect to the chemical potential 
fi related to the charge density of carriers (in the exper- 
iments 0, 0, H, H| , the chemical potential is tunable by a 
gate bias voltage V g ). This phase transition is provided 
by dynamics responsible for creating an excitonic gap A 
in a strong magnetic field: While at small \fi\ the exci- 



*On leave from Bogolyubov Institute for Theoretical Physics. 
03143, Kiev, Ukraine 



tonic gap is generated, there is no gap as becomes 
larger than a critical value |/x c | determined below. In 
fact, if this scenario is correct, the phenomenon discov- 
ered in [8( can be interpreted as the observation of two 
different phases in graphene. 

As will be shown below, one of the predictions of this 
scenario is that the only plateaus in the Hall conductivity 
<j xy are those with v = 0, ±1 and v = ±2fc, k = 1,2,..., 
i.e., the plateaus observed in Ref. 0. The plateau v — 
±1 appears only if both the spin splitting is included 
and the gap A is nonvanishing. Another prediction is 
that the excitonic gap A is much smaller than the gap 
y / 2hv F \eB\/c between the LLL and the n = ±1 LLs. 
In other words, the excitonic gap is produced by weak 
coupling dynamics 0] . This prediction can be checked by 
measuring the critical temperature at which the v = 
and v = ±1 plateaus disappear. We also succeeded in 
reproducing the experimental form of a xy obtained in 

!. 

The fact that a magnetic field is a strong catalyst of 
electron-hole (fermion-antifermion in field theory) pair- 
ing was established long ago (the phenomenon was 
called the magnetic catalysis). The essence of this effect 
is the dimensional reduction D — > D — 2 in the dynamics 
of electron-hole pairing: In a strong magnetic field, this 
pairing is mostly provided by the LLL whose dynamics 
is essentially (D — 2)-dimensional. The dimensional re- 
duction leads to a strong enhancement of the density of 
states. As a result, for D = 2, as in graphene, the pairing 
dynamics in infrared becomes very strong and, for zero 
temperature and zero chemical potential, an excitonic 
gap is generated even at the weakest attractive interac- 
tion between electrons and holes Because of this 
feature, the phenomenon is robust. It was also shown in 
Ref. [10], that at large temperature or/and charge den- 
sity (chemical potential), the excitonic gap disappears. 
As will be discussed below, it also disappears for a large 
impurity scattering rate. 

In graphene, the phenomenon of the magnetic catalysis 
was considered in Refs. [H^Hll m connection with an in- 
terpretation of experiments in highly oriented pyrolytic 
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FIG. 1: Schematic illustration of the spectrum and the Hall 
conductivity in the n = and n = 1 Landau levels for four 
different cases, (a) A = and no Zeeman term, (b) Nonzero 
A and no Zeeman term, (c) A = and the Zeeman term 
is taken into account, (d) Both A and the Zeeman term are 
nonzero. Thickness of the lines represent s the degen eracy 
x4, x2, and xl of the energy states; L — ^/%v 2 F \eB\/ c. 



graphite [14|. In particular, in Ref. the role of 

temperature and chemical potential in this dynamics was 
clarified in detail. Also in that paper, expressions both 
for the diagonal conductivity o xx and the Hall conductiv- 
ity &xy at nonzero gap A were derived and investigated. 
For further studies of this phenomenon in graphene, see 
Refs. [UGl. 

Fig. 1 illustrates the main results of our analysis. It 
shows the spectrum and the Hall conductivity a xy in the 
ii = and n = 1 LLs for four different cases correspond- 
ing to zero (nonzero) gap A and spin splitting. As one 
can see in Fig. Id, when both A and spin splitting be- 
ing nonzero, the plateaus in <j xy observed in Ref. [8| are 
reproduced. Note that the degeneracies of the LLL and 
higher LLs shown in Fig. Id are different. The physics 
underlying Fig. 1 will be discussed in detail in Sees. HT1 
and IIIII below. 

The paper is organized as follows. In Sec. UH we con- 
sider the dynamics of the Hall conductivity in graphene 
when the Zeeman term is ignored (no spin splitting). In 
Sec. IIIII the realistic case, with the spin splitting taken 
into account, is considered. In Sec. IIVI the main results 



of the paper are summarized. In Appendices A, B, and 
C, some useful formulas and relations are derived. 



II. DYNAMICS WITH NO SPIN SPLITTING 

For the description of the dynamics in graphene, we 
will use the same model as in Refs. [12|, HH, the 
so called reduced QED. In such a model, while quasi- 
particles are confined to a 2-dimensional plane, the 
electromagnetic (Coulomb) interaction between them is 
three-dimensional in nature. The low-energy quasi- 
particles excitations in graphene are conveniently de- 
scribed in terms of a four-component Dirac spinor 'F J = 
( , 4'KAa-,ipKBa,il'K'Bcr,^K'Acr) which combines the Bloch 
states with spin a = ±1 on the two different sublattices 
(A, B) of the hexagonal graphene lattice and with mo- 
menta near the two inequivalent points (K, K r ) at the op- 
posite corners of the two-dimensional Brillouin zone. The 
free quasiparticle Hamiltonian can be recast in the "rela- 
tivistic" - like form with the Fermi velocity vf ~ 10 6 m/s 
playing the role of the speed of light: 



H Q = -iv F / d 2 r V a h 1 hV x + 7 2 fiV„) tf e 



(1) 



where ^ a = 4 , J.7° is the Dirac conjugated spinor and 
summation over spin a is understood. In Eq.JT]) 7^ with 
v = 0,1,2 are 4x4 gamma matrices belonging to a 
reducible representation j" = T3 £3 (73, IT2, — vri) where 
the Pauli matrices f , t act in the subspaces of the val- 
ley (K, K') and sublattices (A, B) indices, respectively. 
The matrices satisfy the usual anticommutation relations 
{7^,7"} = 1<r> 9^ = (1,-1,-1), H,v = 0,1,2. The 
covariant derivative V = d+(ie/frc)A includes the vector 
potential in the symmetric gauge A ext = (—By/2, Bx/2) 
corresponding to the external magnetic field applied per- 
pendicular to the plane along the positive z axis. In the 
four-component spinor representation, the Coulomb in- 
teraction has the form 

H int = She J d 2 rd 2 r , ¥CT(r) y)^ (r) _JL_ 

x ¥ ff ,(r)7°vMr), (2) 

where the coupling g = e 2 /e Q hvF and eo is the dielectric 
constant (our convention is e > 0). The total Hamilto- 
nian H to t = Hq + Hi n t possesses U(4) symmetry dis- 
cussed in Appendix [UJ The chemical potential is in- 
troduced through adding the term — /i^y ^ = — /i'l'txj/ 
in H to t (this term preserves the U (4) symmetry) . The 
Zeeman interaction term is included by adding the term 
fiBB^fj cr3^f = ^ f t - 3 ^r ) where now 03 matrix acts on 
spin indices. Here (13 = eh/(2mc) is the Bohr magne- 
ton and we took into account that the Lande factor for 
graphene g^ ~ 2. 

Let us first consider a simpler case with no spin split- 
ting (the Zeeman term is ignored). Then, in App endix 
\K[ utilizing the approach developed in Ref. [17| (and 
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used in [13,[i5|]), we derive the thermodynamic potential 
per unit area in a strong magnetic field B, when the dy- 
namics of the LLL dominate (here the symbol "tilde" in 
the potential and other quantities implies that the spin 
splitting is ignored): 



A/(A,m)- 



Re 



lnT 



1 

7 + i(n + A) 



bL(B) 



/ 2 (A, M ) + 2Tx 



lnT 



2nT 
2nT 



-)-21nr 



2ttT 



,(3) 



where / = y/hc/\eB\ is the magnetic length, L(B) = 
\/hv'p]eB\/c is the Landau scale, T(x) is the Euler 
gamma function, 7 is a LLL impurity scattering rate, 
and the function /(A, fi) is 



1. 



/(A,M) = -Im 



7 + i(n + A) 
2^f 



-(A- -A) 



(4) 

with the digamma function #{x) = ■j-lnT[x). The di- 
mensionless parameter b in Eq. ([3]) reads 



b = 



dk e 



V2j l + fcxo ; 




(5) 



where xo — 0.56\/27rg- The gap equation <9f2/<9A = for 
A takes the form 



A - bL(B)f(A,fi) = 0. 



(6) 



As it is easy to see, for T = fi = 7 = this gap equation 
has only a nontrivial solution, A = bL(B) (the mag- 
netic catalysis). For finite values of these parameters, 
there exists also a trivial solution A = 0, and critical val- 
ues T C ,/^ C ,7 C separate the phases with zero and nonzero 
gaps. The character of the phase transition can be deter- 
mined by studying the thermodynamical potential ([3]) as 
a function of these parameters. Motivated by experimen- 
tal data we are mostly interested in the phase transition 
with respect to the chemical potential /1, because it is eas- 
ily tuned by a gate voltage. The numerical study shows 
that this phase transition can be either a first order or a 
second order one, depending on the values of the scatter- 
ing rate 7 and B. For large enough B (or small enough 
7), it becomes a strong first order phase transition. 

Let us show how the generation of the gap affects the 
form of the Hall conductivity cL. Its expression at 
nonzero A was derived in Ref. [13j |. Here we will use 
a compact expression for a xy , valid for large B and con- 
venient for numerical calculations, obtained recently in 
Ref. El: 



2e 2 sgn(eB) 



'xy 



7t 



2vrT 



ah 

7tr + i(p 



Im yf> 

A) 



7t r + i(fi + A) 1 
2 



2vrT 



2ttT 



(A - -A) ,(7) 



where 7t r is the transport scattering rate (for convenience 
of the readers, the derivation of this expression is pre- 
sented in Appendix IE]) . Note that in graphene, due, for 
example, to a suppression of the backward scattering, 7 tr 
can be smaller than the scattering rate 7 in the thermo- 
dynamic potential ([3]) [l9| . The value of 7t r controls the 
sharpness of the transitions between plateaus in the Hall 
conductivity. 

To illustrate the role of opening the gap, let us consider 
the clean limit case 7 = 7t r = at zero temperature, 
when Eq. takes the form 



2e 5 



sgn/zsgn(e£O0[| M |-A( M ,£)] 



(8) 



When there is no gap, this expression yields the first 
plateau v = ±2 in the half-integer QH effect @, H, 0, H] - 
The appearance of the gap changes the situation: In this 
case, for |/x| < A, the additional plateau a xy = occurs, 
in accordance with the experimental data in Ref. 0. 
Note that the presence of the gap A leads to splitting only 
the LLL. The degeneracy of higher LLs remains the same: 
For these levels , the gap changes o nly the dispersion rela- 
tions, E — ±^j2hvp\neB\/c + A 2 (compare Fig. la and 
Fig. lb). Therefore, besides the plateau v = 0, other 
plateaus are the same as those in the case with A = 0, 
i.e., the filling factors for them are v — ±4(|ra| + 1/2) (see 
Fig. lb). In fact, as will be shown in the next section, in 
this dynamics, the gap disappears for the values of the 
chemical potential corresponding to higher LLs, |n| > 1. 
Then the fact that v = ±4(|n| + 1/2) for these levels 
becomes even more evident. 

Note that lowering the degeneracy of the LLL with gen- 
erating the gap in graphene is connected with the spon- 
taneous breakdown of the initial U (4) symmetry down to 
U(2) a x U(2)b, which is described in Appendix C. The 
point is that the generation of the gap A is connected 
with the order parameter a = = -(0|§*|0) 

[Tfj| . [T2I [l3| and, as shown in Appendix C, it is invari- 
ant only under the U(2) a x U(2)b subgroup of the U(4). 

The following remark is in order. Expression ([7]) does 
not lead to the precise values of the Hall conductivity on 
the plateaus for any finite value of 7t r . This happens be- 
cause localized states are neglected in our consideration. 
With these states taken into account, 7 tr should become 
equal zero between Landau levels (in particular, because 
of that, Eq. ([5]) yields the correct values of the quantized 
Hall conductivity). Thus, strictly speaking, with expres- 
sion iff]), one can describe the Hall conductivity a xy only 
between plateaus, where 7tr controls the sharpness of the 
transitions between the neighbor plateaus. Note how- 
ever that if 7t r is small in comparison with typical energy 
scales in the problem, the plateaus described by Eq. ([7]) 
are rather flat and sharp (see Fig. [3] in Sec. IIIII below V 
As is shown in Sec. IIIII taking 7 tr » 7/3, the behav- 
ior of Hall conductivity described by Eq. ([Jj resembles 
experimental results in Ref. Q. 

Another way to estimate possible values of 7 tr and A 
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FIG. 2: The upper line with triangles: The experimentally 
observed transition points between the lowest plateaus in the 
Hall conductivity from Ref. @]. The lower line: The theoret- 
ical phase diagram in the plane of the magnetic field B and 
the chemical potential for 7 = 18 K and T = 30 mK. First 
and second order phase transitions are denoted by solid and 
dashed lines, respectively. 



is to use the expression 



2e 2 



7tr 



Htt 2 7 2 . + A 2 



(9) 



for the diagonal conductivity derived originally in 
Ref. [l3| and valid at the Dirac point, /i = 0. Although 
for A = it yields u xx = 2e 2 / (?i7r 2 ), which is ir times less 
than the experimentally observed value in Ref. Q , in the 
presence of A, this expression for a xx shows a tendency 
towards an insulating behavior, which is in accordance 
with the data in Ref. d|. 



III. DYNAMICS WITH SPIN SPLITTING 

In the rest of the paper, we will analyze the realistic 
case, with a spin splitting taken into account, whose dy- 
namics is much richer. The thermodynamic potential f2, 
the Hall conductivity a xy and the diagonal conductivity 
a xx are now: 



n(A, M ) = i(f2(A,/i + )+0(A, //_)), 



0+) +a xx (/j,-)), 



(10) 

(11) 
(12) 



where the expressions for f2, a xy and d xx are given in 
Eqs. ((3|), (O, and (|B5j) in Appendix IB] respectively, and 
fi± = /i ± S z , with S z = [IbB being the Zeeman term. 
This term breaks explicitly the U(4) symmetry down to 
U(2) c x U{2) d discussed in Appendix C. 

In order to determine the parameters of the theoretical 
model, we use the experimental data of Ref. [8[ as a guide. 
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FIG. 3: Hall conductivity in the experiment (upper panel) 
and in the theoretical model (lower panel) for magnetic fields 
B = 9 T (circle), 11.5 T (pentagon), 17.5 T (hexagon), 25 T 
(square), 30 T (diamond), 37 T (up triangle), 42 T (down 
triangle), and 45 T (star). The parameters in the model are 
b = 0.04, 7 = 18 K, 7t r = 6 K, and temperature T = 30 mK. 
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FIG. 4: The gap A versus V g for B = 30 T and four different 
values of temperature, T = 30 mK, 5 K, 10 K, and 15 K, from 
top to bottom. 



In the regime when the excitonic gap is sufficiently large 
to remove the degeneracy of the LLL, the model predicts 
that one of the transition points between the plateaus in 
the Hall conductivity corresponds to a phase transition 
in which A — > (see the discussion connected with Fig. [4] 
below) . More precisely, it is the transition point between 
the v = and v — 2 plateaus for B = 11.5 T, and the 
transition point between the v = 1 and v = 2 plateaus for 
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larger values of B (for B — 9 T, with no plateau v = 0, 
the transition point is taken to be zero). This is con- 
trolled by the value of the chemical potential \i (see Fig. 
[5]). In the experiment, the corresponding control param- 
eter is the gate voltage V g . To obtain a relation between 
the two, we compare the theoretical phase diagram on 
the B-fj, plane with the experimental phase diagram on 
the B-Vg plane in Fig. [5] The latter is obtained by a 
simple compilation of the experimentally observed tran- 
sition points between the corresponding plateaus in the 
Hall conductivity. The best fit that was found is: 

,i = 0.5(V g - Va) + 7.0 sgn^ - V )yJ\V g -V \, (13) 

where while /i is measured in kelvins, V g is measured in 
volts, and Vo is the center point in the dependence of the 
Hall conductivity obtained in experiment. Note that the 
values of Vo are different for different values of B and are 
in the interval from 0.8 - 5.8 V [2p|. 

In Fig. [31 we present the experimental data for a xy 
and their description in this model for the values of the 
parameters indicated in the figure. The form of the ex- 
perimental and theoretical curves are quite similar. Note 
that the value of the parameter b — 0.04 corresponds to 
a weak coupling with g = e 2 /eohvF — 0.07. Because of 
the three dimensional nature of the Coulomb interaction, 
it is plausible that the value of g is influenced by a large 
dielectric constant eo of the substrate in the experimen- 
tal device. The main reason of the necessity of a weak 
coupling for the fit is that the lengths of the experimen- 
tal plateaus with v = 0, ±1,±2 imply that the Zeeman 
energy and the gap A are of the same order, and the 
Zeeman energy is only ^10Kat£>^10T. Note that 
although there is no v = plateau in the B = 9 T curve, 
the gap A in this case is also nonzero (although small). 
The reason is that in the presence of a nonzero scattering 
rate 7, the effect of a small gap, A < 7, is unobservable. 

In Fig. |H the gap A versus V g for B — 30 T and four 
different values of temperature is shown. At low T = 30 
mK, a strong first order phase transition with respect to 
V g , i.e. [i, is clearly seen. It corresponds to the transition 
point between the v = ±1 and v = ±2 plateaus. The 
steep descents, which occur approximately at the half of 
the critical value of V g — Vo , are related to the transition 
between the v = and v = ±1 plateaus. The phase 
transition with respect to temperature is a second order 
one with the critical temperature T c ~ 17.2 K. 

Fig. 1 in Introduction summarizes the main results of 
our analysis. Fig. Id clearly shows that an excitonic gap 
A and a large enough Zeeman term (5 Z > 7) together lead 
to the QH plateaus with the filling factors v = 0, ±1 and 
v = ±2k, k = 1, 2,..., i.e., those observed in Ref. [|. Let 
us discuss this point in more detail. It is noticeable that 
while a large enough Zeeman term leads to the plateau 
v = even for A = (see Figs, lc), the v = ±1 plateaus 
appear only if both the Zeeman term is included and the 
gap A is nonvanishing (compare Fig. Id with Figs, lb 
and lc). Therefore the v = ±1 plateaus are the clearest 



signature of the presence of a dynamical excitonic gap. 
There are of course also the plateaus v — ±2 connected 
with the LLL. As was already pointed out in Sec. HjJ the 
degeneracy of higher LLs is removed only by the Zeeman 
term (and not by the gap A). Therefore the filling factors 
of the plateaus with \v\ > 1 are described by v — ±2k, 
k = 1,2,.. ..(see Fig. Id) [2l[. Thus, in this scenario, the 
LLL plays a very special role: While the excitonic gap 
does not reduce the degeneracy of higher LLs, it leads to 
splitting the LLL. This point is at heart of reproducing 
the QH effect data @ in this scenario. 

This picture is intimately connected with the removal 
of the degeneracy in this dynamics. As is shown in Ap- 
pendix C, a nonzero A and the Zeeman term together 
break the initial non-abelian U{A) symmetry down to the 
abelian U(l)\ x £^(1)2 x ( 1 ) 3 x ?7(1)4 one. Since irre- 
ducible representations of abelian symmetries are one- 
dimensional, the t/(4) degeneracy of the LLL is com- 
pletely removed. 

By using Eqs. (TT2"]) and (|B5|) . we also checked that the 
values of A and 7t r utilized in this section yield the diag- 
onal conductivity a xx whose behavior is in a qualitative 
agreement with the data in Ref. Q. 



IV. CONCLUSION 

We believe that the observation of the v = and 
v = ±1 plateaus in the experiment [|[ strongly suggests 
that the existence of a dynamical excitonic gap (or gaps) 
in graphene in a strong magnetic field is a viable pos- 
sibility. In this paper, only a singlet excitonic gap was 
considered. There of course exist other options: For ex- 
ample, the same Coulomb interaction can lead also to a 
nonzero triplet order parameter ("Jer 3 "!'). In that case, 
the states with up and down spins will have different 
gaps, A + and A_, respectively. The same arguments 
as those used above show that the degeneracy of higher 
LLs is lowered only by the Zeeman term and, therefore, 
in that case the filling factors of the the QH plateaus 
will be the same as in the case of the singlet excitonic 
gap. Therefore the modification of Fig. Id will be only 
in replacing A with A± for up and down spins. 

In the present scenario, the weak coupling dynamics 
was utilized. The reason of its necessity is that the 
lengths of the experimental plateaus with v = 0,±1,±2 
imply that the Zeeman energy and the gap A are of the 
same order, and the Zeeman energy is only ~ 10 K at 
B ~ 10 T. Note, however, that this argument is valid 
only for a paramagnetic regime in which there is no large 
enhancement of spin splitting by dynamics in a magnetic 
field. When such a enhancement takes place [22] |. strong 
coupling dynamics might lead to a good fit of the exper- 
imental data in QH effect in graphene. This possibility 
will be considered elsewhere. 

It is instructive to compare the present approach with 
other ones used for the description of the dynamics in 
QH effect in graphene. In Refs. [23|, [24|, the QH ferro- 
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magnetism was considered. The main prediction in Ref. 
[HI is that in this case the QH plateaus with all integer 
values of the filling factor v occur [the critical values of a 
magnetic field B at which plateaus occur are different for 
different v and increase with v\. This prediction is quite 
different from the present one that reflects a difference 
between the dynamics in these two scenarios. As was 
emphasized above, the excitonic gap does not reduce the 
degeneracy of higher LLs and it is unlike the QH ferro- 
magnetism. As a result, there are no odd filling factors 
v = 2k + 1, k > 1; in the scenario with excitonic gaps. 

In Ref. [25| . a scenario with the paramagnetic regime 
was considered (with no large enhancement of spin split- 
ting). Unlike the present scenario, the breakdown of the 
Z7(4) in [25| is not spontaneous but explicit, provided 
by local (on-site) interactions. The main conclusion of 
Ref. [25[ is that Zeeman splitting together with on-site 
interactions can produce QH states at v — 0, ±1 and 
±4 but not at v = ±3 and ±5. Although the values of 
filling factors agree with ours, these two dynamics are 
very different and should lead to very different spectra of 
collective excitations. 

Which of these scenarios is realized in graphene is an 
open issue. It would be interesting to include all possible 
competing orders in the thermodynamic potential and 
find the genuine ground state and the phase diagram in 
graphene [261 ]. 
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APPENDIX A: THERMODYNAMIC POTENTIAL 

In this Appendix we derive expression © for the ther- 
modynamic potential per unit area valid in the strong 
field limit, y/hvp\eB\/c 3> 7, T, /1, A. We will use the for- 
malism of the effective action introduced and developed 



in classical papers [271 128| ] . In those papers, the case of 
the effective action for elementary fields was considered. 
The case of the effective action for local composite fields 
was studied in Ref. [l7| and in our derivation we will 
follow that approach. 

We start from a general definition of the effective action 
in a theory in which the spontaneous symmetry break- 
ing phenomenon is driven by the local composite order 
parameter a = —(^>^) = — (0|WW|0) corresponding to 
the generation of the excitonic gap A [T3, EH . Following 
the conventional way [13, H3, HH , we introduce the gen- 
erating functional W(J) for the Green functions of the 
corresponding composite field through the path integral 

(Al) 

where J(x) is the source for the composite field 
—^(x)^(x) and C qp is the Lagrangian density of quasi- 
particles in the model at hand (note that here x° is the 
time variable t). Then, by definition, the effective ac- 
tion for the field a(x) — — (9 (x)^ (x)) is given by the 
Legendre transform of the generating functional W{J), 



r(cr) = ITT J) - / d 3 xJ(x) cr(x) 



(A2) 



where the external source J(x) on the right-hand side 
is expressed in terms of the field cr{x) by inverting the 
relation 



SW(J) 



a(x). 



(A3) 



5J(x) 

The effective action T(a) in Eq. (|A2[) provides a natural 
framework for describing the low energy dynamics in the 
model at hand. It is common to expand this action in 
powers of space-time derivatives of the field a: 



T(a) = Jd 3 



(A4) 



where V(a) is the effective potential. The ellipsis denote 
higher derivative terms as well as contributions of the 
Nambu-Goldstone bosons. From Eqs. ([A"2|) and (|A3[> . 
we derive the following relation: 



5T(a) 
5a(x) 



-J(x) 



(A5) 



In the limit of a vanishing external source, this equation 
turns into an equation of motion for the composite field 
a(x). In a particular case of constant configurations, the 
equation reads dV/da = 0. 

The thermodynamic potential Q per unit area in Eq. 
© is nothing else but the effective potential V at non- 
zero T and /i. The constant source J plays the role of 
the bare gap (Dirac mass), J = Aq. Then the initial 
relation in the derivation of (following from Eqs. (|A4|) 
and (|A"5|) ) is: 



&7 



(A6) 



7 



At zero T, fi and 7, the gap equation with nonzero Ao 
in a strong field has the form (see Eq.(51) in Ref. [HI)) 



A = Ar 



dto A 



2tt u 2 - A 2 



d 2 k 



■ exp 



fic\k\ 



2eB 



U(k), (A7) 



where 



U(k) 



2tt 



1 



e |k|(l+a|k|) 



a = Am/i 



he 



eohvp y \eB\ 



and the const vq w 0.14 (see Eqs. (46) and (47) in 
At finite T, /1 and 7, the gap equation is written as the 
following sum over Matsubara frequencies uj n — TrT(2n + 
1): 



A = A + 2bL(B)T 



n— — 00 



(cj„+7Sgn(o;„) - ip) 2 + A 2 ' 
The sum over Matsubara frequencies is easily performed, 



r E 



(w„ + 7Sgn(w„) - ifi) 2 + A 2 



= — Im 

2tt 



...,I±^ + i|-(A^-A) 



,(A10) 



where \& is the digamma function. Then the gap equation 
takes the form 



1. 



A = A + bL(B)-lm 



j + i(n + A) 1 



2ttT 



7 + i(M-A) 1 
2vrT 2 



(All) 



Let us express Ao through A from this equation and 
substitute it into Eq. (|A6|1 . Then, taking into account 
the relation 



dil _ dfl dA 
da dA da 1 



we come to the final equation 
dil da 



dA dA 



[A-bL(B)f(A,»)} 



(A12) 



(A13) 



where the function /(A, /1) is given in Eq. (U|). The 
condition dCl/dA = yields gap equation ([6]) in the main 
text. 

Since the field a = — (\f"J r ), we need to evaluate the 
chiral (excitonic) condensate which is given by 



00 . 

(*M>) = -T y, J 



d 2 k 
J^rj 2 



trS(w„, k), 



(A14) 



where the fermion propagator in the LLL is 

S(u n ,k) = 2P_e- k2 ^—— 1 

7o(^n + 7Sgn(w„) - ifjt) - A 

(A15) 

with P_ = (1 -i 7 1 7 2 )/2. 

Calculating the trace and integrating over momenta we 
get 



,- v 2\eB\T ^ 



n— — 00 



(w„ + 7Sgn(w„) - ip) 2 +A 2 ' 
(A16) 

The sum over Matsubara frequencies is evaluated by 
means of Eq. (|A10[) and we obtain 



a = -<**} = ^2 /(A,/i). 



Therefore Eq. (|A13j) can be rewritten as 
0fl _ 1 df(A,n) 



dA irl 2 dA 
Integrating over A we find 

1 



A-bL(B)f(A,n)] 



(A17) 



(A18) 



fi(A,/i) = — ^A/(A, M ) 



bLiB) -f(A^) 



-2TRc 
hlnT 



lnT 



-y + i(p + A) 



2ttT 

7 + iQu-A) 1 
2irT 2 



(A19) 



where the integration constant C(/x) was added on the 
right-hand side. Since C(/z) determines only the overall 
normalization of the potential, we can take it such that 
n(A = 0,/i = 0) = 0. As a result, we arrive at Eq. ([3]) 
in the main text. From Eq. (|A19|) we also find that on 
the solution of the gap equation d^l/dA — 0, the carrier 
density p is given by the expression 



dfl 1 1 
7T = ~P = J2~ lm 

Op TTL Z TT 

7 + f(/x-A) 1 
2ttT 2 



7 + i(/x + A) 
2^f 



(A20) 



APPENDIX B: CALCULATION OF THE 
CONDUCTIVITIES 



In the bare bubble approximation, the expression for 
the diagonal conductivity in the limit of B — > 00 can be 
obtained from Eqs. (3.11), (3.12) in the second paper in 
Ref. 3 (K = 1 in Appendix B) 



e 2 7tr 



du> 



4Tcosh 2 ^ 



(Bl) 



7tr 



7tr 



L7 2 r + (c-A) 2 7 2 +( w + A)2 



8 



[due to reasons pointed out in Sec. [TTJ we use 7t r , and 
not 7, in conductivities a xx and a xy \. The integrals in 
this expression can be evaluated exactly as follows. First 
we write 



derivative with respect to A: 



I = 



Ttr 



cosh 2 £££ 7t 2 r + (w - A) 2 



Re / cfte~ 7tI * 
o 



cosh 2 ^ 



2TRe / dte - *^-*^-^ 



rf.r 



iITU 



cosh a; 



OO OO 

J J cosh x 



(B2) 



daiV e 2 sgn(ei?) 



dA 



4tt 2 T 



Ttr 



da; 



cosh 2 



7tr 



7 2 + (w + A) 2 7t 2 r + (w - A) 2 



Then we use Eq. (|B4|) and find 



(B7) 



The integral over x is evaluated by means of the formula 
(3.982.1) from Rcf. 29]. Then we get 



dtt cos[(fi — A)t] e 
sinh(7rTt) 



-tit, 



(B3) 



dA 



e 2 sgn(e£?) 



2vr 3 T 



Re 



*' 



Ttr + + A) 1 



7tr + »(M - A) 

2ttT 



2ttT 
1 
2 



(B8) 



This integral can be evaluated by differentiating 
Eq. (4.131.3) of Ref. [H[ which yields 



I = 4:7rT Z 



(I 



dttcos[(^- A)t]e~* 7t '- 
sinh(TrTi) 



(B4) 



Thus we obtain 



- _ e 2 7tr R 



2tT 3 T 



VP' 



, f l tI + i(fx + A) 1 
2ttT 2 



*' 



, ^ 7 t r -M(M - A) 1 

2ttT 2 



(B5) 



One can check that for fi = and T — > Eq. (|B5[) reduces 
to Eq. ©. 

Now we derive the expression for the dc Hall conduc- 
tivity. In the limit B — > oo, Eqs. (3.14) and (3.15) in the 
second paper in Ref. [H yield 



Or 



e 2 sgn(e£?) 
4tt 2 T 

7t r (w - A) 



cosh 2 ^ 



+ 



7tr(w + A) 



(B6) 



.7t 2 r + (w - A) 2 7 2 + (w + A) 

cj + A w - A" 
arctan h arctan 



7tr 7tr 

We first consider the terms with arctan functions taking 



Therefore 



e 2 sgn(ei?) 



Im 



7t r + i(/i + A) 



+4- 



Ttr + - A) 
2ttT 



2ttT 
1 
2 



(B9) 



Here we took into account the fact that because of the 
condition a xy (A = oo) = [see Eq. (|B6j) ] and the known 
asymptotics of the vp-function, the integration constant 
equals zero. 

As to the integrals of the two first terms in square 
brackets in Eq. (|B6[) . they are calculated by means of 
the formula 



du> 



cosh 2 7t 2 r + (w - A) 2 



(BIO) 



= Im*' 

7T 



7tr 



-i(M-A) 1 
2ttT 2 



9 



Its derivation is as follows: 

oo 

du> lo — A 

cosh 2 7? r + - A) 2 



Re 



do; 



cosh 2 ^w-A| i 7tr 



where Z4 is the 4x4 Dirac unit matrix and cr Q , with 
a = 0, 1,2,3, are four Pauli matrices connected with spin 
degrees of freedom [a is the 2x2 unit matrix] . The Dirac 
matrices are connected with degrees of freedom reflect- 
ing the band structure of graphene with two inequivalcnt 
Dirac points at the corners of the Brillouin zone. In the 
representation used in the present paper (see Sec. HI]) , 
the Dirac matrices 7 3 and 7 s are: 



Im / dte~ ltA 







dw e i4 (^- A ) 
cosh 2 ^ 



27Tm J dte-^-^-^ 


o 

47Tm / dte-^'-^-^ 



dx 



i2Ttx 



COSh 2 ; 



(Bll) 



dx cos(2Tte) 



cosh 2 : 



AttT 





00 

2 f dtts'm[(n- A)t]e~ t7t ' 



-4ttT^ 



sinh(TrTt) 

OO 

d f dt sm[(fi - A)t] 



G?7t r 7 sinh(7rTi) 




= _^f^+^-A) r 

where, in the very last equality, Eq. (4.131.3) in Ref. [29| 
was used. Combining this contribution with that in Eq. 
(|B9l) . we arrive at Eq. 0. 



7 3 = i 



I 

1 



7 5 = i 



I 
-I 



(C2) 



where / is the 2x2 unit matrix. The order parame- 
ter connected with the generation of the excitonic gap 
A is a = -(**) = -(tfV*) HO, HI HI, where the 
Dirac matrix 7 anticommutes both with 7 s and 7 s and 
commutes with [7 3 , 7 5 ]. The nonzero expectation value 

= {[^ KAa 1p K Aa + tPx-Aa^K'Aa ~ ^ KBa ^KB<r - 

YK'Ba'^K'B ']) ^ S directly related to the electron density 
imbalance between A and B sublattices of the bipartite 
hexagonal lattice of the graphene sheet [H, EH ■ The dy- 
namical generation of the gap leads to the spontaneous 
breakdown of the ?7(4) down to the U(2) a x U(2)t with 
the generators 



o a _ a a 1. 3 5 . 



(C3) 



[note that, as one can see from Eq. \G2\ , [ 7 3 , 7 s ] is diag- 
onal]. As to the Zeeman term, it is connected with the 
spin density operator ^a 3 ^f. Therefore this term explic- 
itly breaks the [7(4) down to the U{2) c x U(2) d with the 
generators 



APPENDIX C: 11(4) SYMMETRY 

The J7(4) symmetry in graphene is discussed for ex- 
ample in Appendix A in Ref. [131 ]. Here we will describe 
the properties of this symmetry used in the main body 
of the paper. 

The 16 generators of the 11(4) are 

a a a a a a a a 1 

— <8>I 4 , ®7 3 , ^-®7 5 , and _ <g, _[ 7 3 , 7 5 ], 

(CI) 



— ^ -2l ( 



>7 



>7 



and 



o a l r 3 5, 



(C4) 

where a' = 0, 3. Eqs. |[C3) and (JUU) imply that the Zee- 
man term and the generation of the gap together break 
the 17(4) down to the XJ(\) X x J7(l) 2 x U(l) 3 x U(l) A with 
the four diagonal generators 



1, 



2 ^[ 7 3 , 7 5 ]; a' = 0,3. (C5) 
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